Surfactant mediated growth of ferromagnetic Mn 5-doped Si 
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We present an investigation of Mn J-doped layers in Si(001) grown by molecular beam epitaxy. 
We discovered that a Pb surfactant has significant effect on the structural and magnetic properties 
of the submonolayer of Mn, which depends on the Si capping layer growth temperature, Tsi, and 
the Mn coverage, 0Mn- The results presented in this paper identify three regions in the growth- 
phase-diagram characterized by distinct magnetic behaviors and crystal structures. In one region, 
X-ray absorption fine structure (XAFS) and transmission electron microscopy (TEM) experiments 
indicate that MnSi nanocrystallites form with B2-like crystal type structures. At the optimal growth 
conditions, Tsi = 200 °C and 9mh — 0.26 monolayer, a ferromagnetic phase develops with a Curie 
temperature Tc > 400 K and a saturation moment m 3at — 1.56 /is/Mn, whereas Tc drops to zero 
for a control sample prepared without Pb. For Tsi > 200 °C, MnSi-B20 type precipitates form with 
a Tc ~ 170 K. Rutherford backscattering spectroscopy shows that the increase in the remanent 
magnetization in these two phases is possibly correlated with an increase in the Mn substitutional 
fraction, which suggests that a Sii-xMnj, dilute magnetic semiconductor may be forming in the 
matrix between the precipitates. Density functional calculations show that Pb changes the pathway 
by which the Mn atoms access the Si substitutional sites, Mnsi. While the Pb increases the formation 
energy of Mnsi at the Si surface, it enables substitutional incorporation by lowering the formation 
energy of Si vacancies by 0.92 eV. 

PACS numbers: 75.70.Ak, 61.05.cj, 75.50.Pp, 85.75.-d 



I. INTRODUCTION 

Dilute magnetic semiconductors (DMS) offer control 
over the carriers responsible for magnetic order, provid- 
ing unique opportunities for modern electronics* 1 - Given 
the technological prevalence of Si, it is important to un- 
derstand whether it is possible to create a Si-based DMS. 
Mean field calculations predict relatively high Curie tem- 
peratures, Tc, for a dilute Sii-^Mn^ alloy^ and density 
functional theory (DFT) calculations show that substi- 
tutional Mn in Si will produce a ferromagnetic semicon- 
ductor with a magnetic moment just below 3 hb /Mn. 3 
One of the challenges that Si presents for DMS is the 
high growth temperatures that are required for epitaxy^ 

Early attempts to create a homogeneous alloy by 
molecular beam epitaxy (MBE) demonstrate that the Mn 
phase segregates into nanocolumns£ In Mn-implanted Si, 
post implantation annealing produces precipitates from 
the MnSii.7 family of phascsjS However, the nanocrys- 
talline phases that result are interesting in their own 
right. Sputtered Mno.05Sio.95 alloys^ and Mn-implanted 
gj8j9 rC p 0rTj Curie temperatures in excess of 400 K, in 
contrast to all known bulk Mn-Si phases that all have 
magnetic ordering temperatures below 47 K. The origin 
of this high-T C ' phase is not clear (see discussion in Ref. 
[Tol and references therein). The magnetic properties that 
are sensitive to the preparation conditions ^have been 
suggested to be due to the presence of Mn impurities 



in the MnSii.7 phased A recent proposal attributes the 
high Tc to the interactions between the large magnetic 
moments of the impurities inside the precipitates The 
Si matrix must play a role in coupling these particles, al- 
though the exact mechanism is still not well understood. 
A second interpretation suggests that carrier-mediated 
ferromagnetism due to Mn impurities in the Si matrix is 
responsible for the high Tc<£ 

Another way to create a DMS is by the deposition of 
submonolayers of magnetic impurities in a semiconduc- 
tor host, referred to as (5-doping, as first demonstrated 
in (Ga,Mn)As and (Ga,Mn)Sb^ii DFT calculations 
show that a <5-doped layer of substitutional Mn in Si is 
half-metallicjiii However, scanning tunneling microscopy 
studies show that interstitial sites arc energetically more 
favorable^ which drives Mn to the subsurface where it 
forms wiresi 17 i 18 Interstitial sites are also interesting can- 
didates for creating ferromagnetic order since DFT cal- 
culations show that a 1/4 (5-doped layer of interstitial Mn 
in Si is half-metallicJ^ A pure interstitial layer is likely 
difficult to achieve, due to the tendency of the Mn to 
form clusters ; 20 ' 21 although a mixed (5-doped layer may 
also yield interesting magnetic properties^ 

In this paper, we investigate the influence of a Pb sur- 
factant on the local structure and magnetism of (5-doped 
Mn-layers on Si(001). In the case of Si, <5-doping of- 
fers the advantage of independent control of the Mn and 
Si growth temperatures, where the Mn can be grown at 
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FIG. 1. (color online) Growth phase diagram based on Mn 
coverage 0Mn and Si cap growth temperature Tsi. The three 
control samples without Pb are shown in red circles. We iden- 
tify three regions: A, B, and C in the growth diagram based 
on the magnetic and structural properties of the samples dis- 
cussed in the text. The dashed line bounds a mixed-phase 
in region-C. We define 1 ML to correspond to the density of 
atoms in the Si(OOl) plane (6.78 x 10 19 atoms/m 2 ). 



0.86 ML and Si-cap growth temperatures 150°C < T Si < 
300°C. The three regions of the growth phase diagram 
are presented in Fig. [TJ labeled A - C, and the sample 
descriptions are listed in Table [TJ Superconducting quan- 
tum interference device (SQUID) magnetometry mea- 
surements show three distinct regions of magnetic be- 
haviors. In region-A, we create a material with Tq > 
400 K that results from two cooperative phases, where 
one is a nanocrystallinc B2-likc MnSi precipitate and the 
other is possibly from Mn occupying Si substitutional 
sites (Mnsi). In region-B, ferromagnetism is absent. In 
region-C, large nano-disks with a B20-like MnSi struc- 
ture precipitate to form a magnetic phase with an order- 
ing temperature Tq ~ 170 K. Density functional theory 
(DFT) calculations are presented to give insight into the 
role of the surfactant. Lead is found to decrease the for- 
mation energy of Si vacancies, and thereby create a mech- 
anism for increasing the concentration of substitutional 
Mn. 



II. METHODS 



cold temperatures and the Si cap at higher temperatures. 
In this way we attempt to avoid the large Mn precipi- 
tates that form during the annealing of submonolayers of 
Mn deposited on Si(001).— By adding two atomic lay- 
ers of Pb that surface segregate during growth, the sur- 
face diffusion length of Si increases^ which enables the 
crystalline growth of the capping layer at a lower sub- 
strate temperature than is otherwise possible on a clean 
Si(001) surface.— Lead also has the advantage of being 
a group IV clement that does not act as a dopant in Si 
and has been successfully used to incorporate As in Si at 
low temperatures The presence of the surfactant also 
modifies the formation energies of the dopants, as found 
in GaP;2Zr— which could give rise to other important op- 
portunities. 

Our previous x-ray absorption fine structure (XAFS) 
results, however, show that the Mn <5-doped layers 
also have the strong tendency to form precipitates^ 
Previously^ we showed that the structure of ultrathin 
films of Mn on Si is highly sensitive to Mn concentration 
and to substrate temperature during both the Mn growth 
and the growth of the Si capping layer. Our XAFS and 
magnetometry experiments showed that the growth tem- 
perature of the Si capping layer influences the local coor- 
dination of Mn and leads to MnSi precipitates with B2- 
and B20-like crystal structures with completely different 
magnetic properties, and a T c varying between 10 K to 
above 250 K. The B2-MnSi phase was predicted by DFT 
calculations to have a high Tc*^ 

With the complementary data obtained from XAFS, 
Rutherford backscattering spectroscopy (RBS) and 
transmission electron microscopy (TEM) measurements, 
we investigated the structural evolution of the <5-dopcd 
layers in a range of Mn coverages 0.26 ML < 6>Mn < 



A. Growth 

We used MBE to grow Mn <5-doped films on low (1-20 
O cm) and high resistivity (600-1200 ft cm) boron-doped 
Si(001) substrates. The Si wafers were prepared following 
the method of Ref. [U For all samples, a 20 nm Si buffer 
layer was deposited onto the (2x1) surface- reconstructed 
Si(001) prior to thin film growth. We deposited 2 ML 
Pb, measured by reflection high-energy electron diffrac- 
tion (RHEED) intensity oscillations, at room tempera- 
ture. By the completion of the second layer, a (2 x 1) 
reconstruction of the Pb covered surface was observed 
againi 33 i 34 



TABLE I. Samples and their corresponding 8me, Tsi and 
region in growth-diagram. We refer to the samples in the 
text as listed in the first column. 



Sample 


0Mn (ML) 


T si (°C) 


region 


ai 


0.26 


200 


A 


0,2 


0.32 


200 


A (No Pb) 


bi 


0.59 


200 


B 


b 2 


0.68 


200 


B 


bz 


0.7 


200 


B (No Pb) 


Cl 


0.26 


250 


C 


C2 


0.41 


250 


C 


C3 


0.49 


250 


C (No Pb) 


C4 


0.86 


250 


C 


C5 


0.33 


300 


C 



The substrate temperature was lowered to T < °C 
before the deposition of a submonolayer of Mn. The Mn- 
flux rate was monitored during growth with a Bayard- 
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Alport ionization gauge placed in the Mn beam, and 
the amount of deposited Mn was subsequently confirmed 
with RBS. Ex-situ x-ray photoelectron spectroscopy mea- 
surements confirm that Pb segregates to the surface dur- 
ing the growth of the 12.4 nm-thick Si capping layer 
that was deposited onto a substrate held at a temper- 
ature Tgi- The RHEED diffraction pattern showed that 
the cap remained crystalline throughout the growth, al- 
though TEM shows stacking faults in the capping layer 
in samples in rcgion-B and region- C. We also prepared 
control samples without a Pb-layer in order to determine 
how the Pb affects the Mn incorporation and the mag- 
netic properties. 



Energy (MeV) 
1.0 1.2 1.4 1.6 



2.0 



RBS 



We determined the Mn-dopant concentration and the 
fraction of Mn atoms that occupy Si substitutional sites 
(^sub) from RBS. The 2 MeV He ions were directed along 
the [100] direction of the sample. A wide angle detector 
(WAD) measured the ion yield backscattered through an 
angle of « 170°, and a glancing angle detector (GAD) 
collected ions with grazing incidence. Since the WAD 
measurements are less prone to small changes in geome- 
try, we used those values to calculate the amount of Mn 
and x su b, whereas the values obtained from GAD were 
employed to determine the depth of the Mn layer. We 
estimated x S ub from the difference in the Mn signal mea- 
sured in random and channeling orientations. Usually, 
the Mn concentrations determined from random orien- 
tation spectra are more reliable than the Mn amounts 
determined from channeling orientation because they are 
compared to the Si random height, whereas the channeled 
spectra depend on less reliable charge normalization to 
extract the Mn concentration. Therefore we checked the 
values of x su b by using the signal from Pb as a reference, 
assuming that the Pb peak does not allow channeling. 
As a further check on the reliability of charge integra- 
tion, and specifically for control samples without Pb, we 
verified that the minimum ion yield measured from the 
spectrum height from ions scattered deep in the Si was 
within the expected range of 4-5% for a well channeled 
spectrum. Furthermore, since in a channeling experiment 
the incident ions tend to travel in the center of the chan- 
nel, the probability of backscattering from interstitial Mn 
is enhanced, this is referred to as flux peaking. Because 
of flux peaking, the channeled signal overestimates the 
amount of interstitial Mn and, consequently, the calcu- 
lated substitutional Mn fractions should be interpreted 
as lower thresholds values^ Characteristic RBS spectra 
measured in random and channeled orientations and the 
variation of x su b with #Mn are shown in Fig. [2] 




250 300 350 400 450 500 550 600 
Channel 




FIG. 2. (color online) (a) RBS spectra of one of the samples 
collected in random and channeled orientations. The Mn sub- 
stitutional fraction x su h is calculated by using the difference 
between the areas of Mn peaks in both orientations, (b) Vari- 
ation of a; su b with Mn coverage #Mn for samples in all three 
regions of the growth diagram. Open triangle (circle) refers 
to control sample of region-A (region-C) grown without Pb 
and the dashed line is a guide to the eye. 



C. XAFS and TEM 

We probed the local coordination of Mn by XAFS in 
fluorescence mode. Polarization-dependent XAFS spec- 
tra were collected at the Mn K-edge with the samples 
mounted on a spinner turning near 30 - 40 Hz to minimise 
contributions from Bragg peaks. X-rays were incident on 
the sample at an angle between 5 and 7° above the plane 
of the wafer (83 - 85° off the surface/substrate normal). 
Hence, for out-of-planc (oop) polarization measurements, 
the electric field vector of the linearly polarized X-rays 
was then oriented between 5 - 7° off the [001] substrate 
normal direction of the spinning wafer. In-plane (ip) 
polarization measurements were averaged over all in- 
plane orientations perpendicular to [001]. Given the high 
symmetry of the [001] direction in Si, no orientation- 
dependence was expected in the in-plane XAFS signal. 
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The XAFS experiments and analysis were conducted as 
outlined in Ref. [HI using the facilities of the PNC/XSD 
ID beamline 3 -^ at the Advanced Photon Source, sector 
20. 




1 2 3 1 2 3 4 5 



R (A) R (A) 

FIG. 3. (color online) (a) In-plane XAFS data (insets) and 
corresponding fit in i?-space for samples a\ and C4. Profile 
functions of Mn-Si and Mn-Mn scattering paths are shown for 
each sample and the fitting range is within R = 1.40 -2.82 A. 
In (b) and (c) we compare the \FT[k 2 x(^)]| of samples ai and 
c 4 with the calculated \FT[k 2 X (k)}\ of MnSi-B2, MnSi-B20 
and Mn4Si7 structures. We used the same window function 
and similar fc-range for the FT of the calculated structures 
and scaled their \FT[k 2 xW]| to match either the Si (MnSi- 
B2 and Mn 4 Si 7 ) or the Mn (MnSi-B20) peak of the measured 
spectra. 

We obtained the coordination parameters from fit- 
ting the magnitude of the Fourier transform (FT) 
of the fc 2 -weighted XAFS interference function x(^)> 
\FT[k 2 x{k)]\, in i?-space. We show two XAFS xO) in 
Fig. [3] that are representative of our XAFS data. As can 
be seen from the insets of Fig. EJa), there is an increase 
in the noise level of k 2 x(&) for k > 10.0 A -1 therefore, 
we limited the FT to within the closest zero-crossing of 
the fc-rangc = 2.3 — 10.0 A -1 . We included two scatter- 
ing paths that corresponded to the first Mn-Si shell and 
the second Mn-Mn shell, with the exception of sample 
61-ip, which had two separate Mn-Si shells. To further 
confirm the composition of the two shells, we applied the 
equations derived from beats theory in Refs. [33 and EjH 
to calculate the difference in radii between the first two 
shells, Ai? = R 2 — Ri. The values we found were con- 
sistent with the Ai? that we obtained from our best fit 
values of i?i and i?2. 

In Table [U we report the XAFS parameters of the 6- 
doped samples. For comparison, we also include the pa- 
rameters of a 5 nm MnSi/Si(lll) reference sample with 
a B20 crystal structure, as well as the parameters ex- 
pected for Mn in a Si substitutional site and for a Mn in 



a tctragonally distorted B2 (CsCl) structure. The am- 
plitude reduction factor Sq = 0.70 and the energy shift 
AEq = 5.9 eV were fixed to the values obtained from 
the reference MnSi-B20 sample. Given the weak Mn sec- 
ond shell and the strong correlation between N 2 and the 
mean square relative displacement a\ for the measured 
samples, we could not determine a value of N 2 - Instead, 
we refer to the area under the Mn-Mn profile function, 
■A-Mri—Mm to represent the relative combined change in 
N 2 and of. 

We further characterized the structure of the samples 
with plan view and cross-sectional TEM, imaged with 
300 kV Philips CM30 TEM. Samples were prepared by 
low-angle mechanical polishing and subsequently cleaned 
with Ar-ion milling^ We also used energy dispersive x- 
ray spectroscopy (EDX) to probe the chemical composi- 
tion of the precipitates that formed in certain samples. 



D. Magnetometry 

We analyzed the magnetic properties of the samples by 
measuring their moment as a function of field (m-H) and 
remanent moment as a function of temperature (m r (T)) 
with a Quantum Design MPMS XL-5 SQUID magne- 
tometer. For all magnetic measurements, 3-5 pieces of 
each sample were mounted inside a sample holder con- 
sisting of a clear plastic straw. We did not find a differ- 
ence in the magnetization when measured with the field 
applied parallel or perpendicular to the [110] orientation 
of the films, which is likely due to disorder in the Mn 
layer. To remove the paramagnetic and diamagnetic re- 
sponses of the low-resistivity Si(001), where the boron 
concentration is of the order of 10 15 — 10 16 atoms/cm 3 , 
we measured the magnetic response of a control sam- 
ple that consisted of a 100 nm Si buffer layer grown on 
Si(001) and subtracted it from that of the <5-doped sam- 
ples. We could not detect a paramagnetic signal from a 
similar control sample grown on high-resistivity Si(001), 
and therefore only the diamagnetic response of the Si was 
subtracted from the samples grown on the high-resistivity 
substrates. Furthermore, we measured the magnetic re- 
sponse of a sample prepared with 2 ML Pb under the 
same conditions of the <5-doped films but without Mn. 
We found that the saturation and remanent moments 
of this Mn-free control sample are zero within the error 
of the measurement, which confirms that the remanent 
magnetization in the 5-dopcd films arises from Mn. 

We calculated the Mn saturation moment from the m- 
H loops measured at T = 2 K. For m r (T) measure- 
ments, we saturated the magnetization in a 5 T applied 
field, and then to ensure that the field is set to zero, we 
used the magnet reset option of the SQUID magnetome- 
ter to discharge the magnetic flux trapped in the coils. 
Measurements of m r (T) were then made while increas- 
ing the temperature from 2 K. We estimated the Curie 
temperature Tq from the point where m r (T) drops to 
zero. 
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TABLE II. Coordination numbers (N) and distances (R) obtained from the fit of the in-plane (ip) and out-of-plane (oop) XAFS 
data of the samples and the ip data of a 5 nm MnSi-B20 film. Estimated parameters of a tetragonally distorted t-MnSi-B2 
structure and of Mn in Si substitutional site, Mnsi, are also included. Awn-Mn is the area under the Mn-Mn profile function 
relative to that of Mn-Si. The percent residual (R%) and uncertainty of each fitting parameter (in parentheses) are also 
included. The energy shift AEo = 5.9 ± 0.6 eV and S% = 0.70 ± 0.02 were obtained from the reference MnSi-B20 sample and 
treated as fixed parameters in the 5-doped samples. 



Samples 


AT ( Qa\ 

iVi \b%) 


Hi (A) 


a x (A ) 


AT ( 1\/Trr^\ 


rl2 (A) 


°2 (A ) 


A (Of \ 
AMn-Mn (/o) 


r? (Of\ 
ti \/o) 


ai-ip 


6.7(7) 


1.6 1(1) 


U.UU9(2 ) 




z. ) 




ICQ 


J..U4 


ai-oop 


7(1) 


2.378(5) 


0.010(2) 




2.75(3) 




15.0 


1.14 


&i-ip 


2(1) 


2.22(3) 


0.006 


7(l)(Si) 


2.42(2)(Si) 






2.61 


fei-oop 


6.1(8) 


2.38(1) 


0.008(2) 




2.80(1) 




21.4 


0.91 


&a-ip (No Pb) 


7.7(4) 


2.365(5) 


0.0096(8) 




2.79(1) 




9.2 


0.47 


& 3 -oop (No Pb) 


7.0(6) 


2.367(4) 


0.008(1) 




2.82(2) 




13.2 


0.82 


c 4 -ip 


6.3(6) 


2.38(1) 


0.010(2) 




2.817(5) 




41.0 


1.11 


C4-OOp 


6.3(6) 


2.38(1) 


0.012(1) 




2.829(4) 




44.9 


0.62 


MnSi-B20 


6.5 


2.39(1) 


0.010(4) 


5.9 


2.788(6) 


0.005 




1.88 


t-MnSi-B2 


8 


2.370 




6 


2.790 








MnsC 


4 


2.351 















a We consider 1 substitutional Mn atom in the Si lattice since 6>Mn < 1 ML. In this case, the second shell corresponds to Mn-Si with 
R2 = 3.84 A, which is outside our fitting range and therefore not included in the table. 



III. MAGNETIC AND STRUCTURAL 
PROPERTIES 

SQUID results demonstrate distinct magnetic behav- 
iors in each of the three regions of the growth-phase di- 
agram shown in Fig. [TJ Region- A has a Tc well above 
room temperature, as indicated by the non-zero rema- 
nent moment measured up to T = 400 K (Fig.|4ja)). The 
existence of magnetic order is confirmed by the presence 
of hysteresis shown in Fig. EJa). Sample a\ was found 
to have the optimal growth conditions that produced the 
largest m r (Fig. [4] (a)) and the largest saturation mo- 
ment, m sat = 1.56 ^b/Mii (/is thereafter) measured at 
T = 2 K (Fig.Efd)). This sample has the lowest Mn con- 
centration, #Mn = 0.26 ML, and a growth temperature 
that is closest to the phase boundary between regions-A 
and C. The m r decreases as Tg; drops from this optimal 
Tg; = 200 °C (Fig. Bid)), while the m sat measured at 
T = 2 K drops quickly with increasing 6*Mn (Fig. [Hd)). 

Fits to the XAFS spectra resolve the first two shells 
surrounding the Mn atoms that correspond to the Si 
nearest neighbors and Mn next-nearest neighbors (see 
Table HJ). The Mn-Si bond length, i?i = 2.374 A and 
the coordination number Ni = 6.7 ± 0.7 are larger than 
expected for substitutional Mn (see Table |LT|) . The large 
coordination indicates the formation of precipitates, al- 
though these phases could not be detected in our TEM 
experiments. A qualitative comparison between XAFS 
measurements of sample a\ and calculated spectra shows 
that the local Mn environment does not correspond to 
either MnSi-B20, MnSii.7 or MnsSi3 crystal structures 
that arc commonly produced from deposition of Mn 
on Si(001) i 40 ' 41 Furthermore, the Mn-Mn bond length 



i? 2 = 2.75 A is significantly contracted relative to that 
measured for MnSii.7 precipitates in Ret. [I2T 

According to DFT predictions, tetragonally distorted 
MnSi-B2 is the most energetically favorable structure 
that forms on Si(001), since MnSi-B2 phase has a rea- 
sonably good lattice match to Si(001)42 We reported the 
existence of a B2-like phase in our previous work of MnSi 
films on Si(001)4i A comparison between the fitting pa- 
rameters in Table [TTI with those expected for a B2 show 
that this structure is the closest match to sample 01. 
The expected 8-fold coordinated crystal structure can be 
constructed from a Si lattice by placing Si and Mn into 
interstitial sites. Missing intcrstitials would account for 
some of the discrepancy between the measured and the 
expected coordination number for a B2 structure. RBS, 
however, shows that only 65% of the Mn in sample 01 is 
in interstitial sites, which indicates that Mn in this sam- 
ple has at least two distinct environments. Although the 
XAFS is not inconsistent with the co-existence of both a 
B2 phase and Mn in Si substitutional sites (Mnsi) such a 
fit would introduce too many fitting parameters to per- 
mit a quantitative analysis. 

The SQUID data is also suggestive of two components: 
FigJSJa) reveals a wasp-waisted hysteresis loop that can 
be constructed from the superposition of a ferromagnetic 
and a superparamagnetic hysteresis loop, as shown in the 
case of MnSii.7 nanoparticlcs in Si4£ The ferromagnetic 
hysteresis loop seen at low magnetic field in the inset 
of Fig. 02a) has a weaker temperature dependence than 
the s-shaped portion of the m — H curve visible at high 
fields, which suggests that these components have differ- 
ent magnetic behaviors. 

The samples in Region B of the growth-phasc-diagram 



6 



(a) 



-A-0.33 ML, T s =150°C 
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FIG. 4. (color online) (a) The non-zero m r (T) indicate a 
ferromagnetic phase with Tc > 400 K for samples in region 
A. The control sample of this region (filled triangles) has zero 
remanent moment, (b) m r (T) show no ferromagnetic ordering 
for samples grown with Pb in region B but a Tc > 400 K 
for the control sample, (c) A magnetic phase transition is 
observed in region-C with Tc — 170 K for all samples with 
Pb. There is no remanent moment for the control sample of 
this region, (d) Remanent moment measured at T = 250 K 
as a function of #Mn for all Tsi- Dashed lines connect the 
samples grown with same Tsi. 



have Tc's near zero (Fig. |4jb)) and the m sa t values are 
small (Fig. E^d)). However, there is a noticeable change 
in the Mn environment when 9mh increases to 0.59 ML 
at Tsi = 200 °C. In particular, we note the distinctive in- 
plane structure of sample b\ which has two Mn-Si shells 
as opposed to Mn-Si and Mn-Mn shells for all other sam- 
ples. The short Mn-Si distance R\ = 2.22 A forms only 
in plane. The TEM images of Fig.(6^a) show that sample 
62 contains precipitates with an estimated lateral size of 
up to 15 nm. Curiously, samples from this region still 
show a considerably large x su h (Fig. [5]). The small m sat 
and low Tc indicate that these sites are not magnetically 
active. 

In region-C of Fig. [TJ the m r has a concave shape typ- 
ical of ferromagnetism that drops either to zero or to 
a small constant value at a temperature of 170 K. This 
transition temperature is unaffected by a variation in 0Mn 
or an increase in T si to 300 °C (Fig. [H c ))j which reflects 
a higher stability relative to the phase in region-A. It is 
likely that this magnetic behavior arises from the large 
precipitates that form in this region and can be seen in 
the cross-section TEM image of sample C4 (Fig. EJc)), 
given that the onset of this new low temperature feature 
coincides with their appearance. The nano-disk shaped 
precipitates (Fig. Etd)) distribute themselves at the in- 
terface with a range of sizes that vary between 10-200 nm 
in diameter. 

Energy-dispersive x-ray spectroscopy (EDX) showed that 



FIG. 5. (color online) m-H loops measured at low and high 
temperatures as indicated. The insets show that there is a 
hysteresis at T = 300 K for samples (a) a\ and (b) c\ but 
it vanishes for sample C4 (c). (d) Variation of Mn satura- 
tion moment m sa t calculated at T = 2 K as a function of 
#Mn- Open symbols refer to samples grown on low resistiv- 
ity Si(001) where the B paramagnetic contribution has been 
subtracted. 



the precipitates contain Mn, but it was not possible to 
determine their phase from the weak reflections in the 
selected area diffraction pattern (SADP) measurements. 
A comparison between the XAFS function of sample C4 
and that for the MnSi-B20 structure shows similar fea- 
tures up to R= 5 A (Fig. Etc)). Although the first Mn-Si 
shell coordination number and radius are within the error 
limits of those of MnSi-B20, the ip and oop bond lengths 
of the second shell appear to be expanded relative to 
the i?2 of a B20 structure (see Table Hi]) . Furthermore, 
the transition at T = 170 K is not consistent with the 
T c < 44 K of ordered MnSi-B20 thin films 44 We expect 
that the low 0Mn and low Tsi create a defected MnSi-B20 
structure. Such structures have been observed in similar 
Mn-doped Si systems. 31 ' 45 

The presence of a second phase is evidenced by the 
non-zero m r above 170 K in the lower #Mn samples in 
region-C. The m r of sample c\ above 170 K matches 
the temperature dependence of the remanent moment of 
region-A samples. Furthermore, the values of the high- 
temperature m r for the region-C samples are comparable 
to those of region-A (see Fig. BJd)) with the exception 
of the optimal sample a\. This suggests that while the 
higher Tgi drives the Mn into larger precipitates with 
Tc ~ 170 K, some of the Mn forms the same structure 
found in region-A. Like the phase in region-A, this second 
phase is sensitive to 6*Mn- Figure |4jc) shows that this sec- 
ond phase is absent from sample a (#Mn = 0.86 ML). It 
is also sensitive to growth temperature: m r (T) = above 
T = 170 K for the T si = 300 °C sample in Fig. He). 

The contrasting magnetic behavior of samples in the 



50 nm 



FIG. 6. (a) The bright field plan-view TEM image of sample 
ai hardly shows any precipitates whereas (b) Mn-rich pre- 
cipitates are clearly observed in sample 62. (c) Cross-section 
TEM image of sample C4 shows Mn-rich precipitates (dark 
features indicated by arrow) distributed at the interface be- 
tween the substrate (1) and the Si-cap (2). The epoxy layer 
(3) that glues the samples to a second Si wafer (4) is also 
visible in the image. The precipitates in this sample have 
nano-disk shape as seen in (d) 



three regions of the growth-phase diagram raises the 
question of the mechanism of the magnetic order. For 
regions- A and C, the m r and the m sa t are found to 
decrease with increasing #Mn, as shown in Figs. |4] (d) 
and[SJd). As the concentration of Mn increases, so does 
the probability of the formation of precipitates. It is in- 
teresting that x su b also follows this trend (Fig. [2|) for 
regions- A and C. Given the nanoscale size of the precip- 
itates, however, it is important to understand the role 
that the Si matrix must play in mediating the ferromag- 
netic order. One possible mechanism is the formation of a 
DMS in the regions between the precipitates. RBS mea- 
surements show that the remanent magnetization is pos- 
sibly correlated with the fraction of substitutional Mn, 
Xsub- In Fig. [71 the m r increases approximately linearly 
with x su h for samples in regions-A and C. RBS measured 
%sub = for sample c\ , although flux peaking may be un- 
derestimating this value, and the uncertainty in the mea- 
surement places the value 1.5 standard deviations from 
the line of best fit. Since MnSi-B2 samples are known to 
be ferromagnetic with a high T C r^ a Sii-^Mn^ DMS 
provides a plausible explanation for how the precipitates 
in region- A are correlated in order to yield a high Tq- 
A similar mechanism of DMS mediated interaction ex- 
plains the coupling between the Mn^Gci-j nanocolumns 
in Ref. l46l and was shown to compete with an antiferro- 
magnetic dipolar coupling between the nanocolumns. 

In region-C, however, the MnSi-B20 precipitates were 
not expected to have a Tc as large as 170 K due to the 
small Tq < 44 K found in MnSi thin films4£ It is pos- 
sible that the model put forward in Ref. [l^ to explain 



FIG. 7. (color online) Variation of remanent moment mea- 
sured at T — 250 K with the fraction of substitutional Mn, 
a; su b. A possible correlation exists for samples of region- A 
(blue triangles) and region-C (red circles), as shown by the 
black line. 

the high Tq in nanocrystalline MnSii.7 precipitates in 
Si also explains the behavior of the region-C samples. 
Like MnSii.7, MnSi-B20 is a weak itinerant magnet. The 
defects suggested by the XAFS data may create large 
localized magnetic moments that are coupled by spin- 
fluctuations in the MnSi, as proposed by Men'shov et 
alJ£ A DMS may be the second phase that is mediating 
the interaction between these large moment precipitates. 
The explanation could also be related to the interfacial 
Mn, which is shown theoretically and experimentally to 
have an enhanced magnetic momentj 42 ' 44 . It is impor- 
tant to note that, unlike a similar structure of Mn doped 
Ge where the Mn acts as a surfactant and floats to sub- 
surface sites in the Ge cap 4 ^, our measurements show 
that Mn in the (5-doped Si layers remains localized at the 
interface. For instance, we could not detect Mn-rich in- 
clusions in the Si cap from the cross-sectional TEM image 
of Fig-EJJc), but we only observed the MnSi nano-disks at 
the interface. Also, the depth of the backscattering Mn 
atoms measured by RBS GAD detector was consistent 
with the thickness of the Si capping layer (<~ 12.4 nm), 
which indicates that Mn is buried at the interface. 

While Fig. [7j shows a possible correlation between m r 
and a; su b, we observed no correlation between m sat and 
a; su b that would allow us to extract the magnetic mo- 
ments of the individual magnetic phases in our composite 
samples. A difficulty with the interpretation of the RBS 
data is that we cannot distinguish the substitutional-like 
Mn in the precipitates from the Mns;. 

Surprisingly, the region-B samples were found to have 
the highest x su h (Fig. EJ^b)). The absence of ferromag- 
netism in the region-B samples could be due to a num- 
ber of possibilities. One possibility is the presence of 
an antiferromagnetic component. DFT calculations of 
Mn interstitial-substitutional complexes, however, pre- 
dict ferromagnetic behaviour^ More recent calculations 
of (5-doped layers of mixed substitutional and interstitial 
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Mn finds a complex non-collincar magnetic order^ A 
second possibility is that the measured x su h of the sam- 
ples of region-B arises mostly from the precipitates and 
a lack of Mngj weakens the carrier mediated exchange 
interaction, which becomes dominated by a possible an- 
tiferromagnetic dipolar interaction between the precipi- 
tates. Another example of a low moment Mn-Si phase 
is amorphous Sii-^Mn^ alloys^ where it has been sug- 
gested that the Mn moment is lost in an itinerant band 
created by Anderson localization. 



IV. THE INFLUENCE OF THE PB 
SURFACTANT 

The control samples grown without Pb demonstrate 
the dramatic influence of the Pb surfactant on the mag- 
netic properties of the Mn (5-doped layers. Without Pb, 
the samples in regions-A and C have no long-range mag- 
netic order, as shown by the zero remanent magnetization 
of these samples in Fig. IUa) and (c), although the satu- 
ration magnetization for the control sample in rcgion-A 
is 0.78 ug. Surprisingly, the situation is reversed in re- 
gion B. The sample without Pb has a Tq > 400 K and 
a large remanent magnetization (Fig. H(b)). This behav- 
ior is explained in our previous work where a MnSi-B2 
structure forms from the deposition of Mn on Si(001).— 
It is, however, worth pointing out that this B2 phase 
grown without Pb only forms at the higher coverages rep- 
resented by region-B, and RBS indicates that the x su h is 
close to zero for this sample. Like the other samples, we 
find that the x su b measured by RBS is not a good in- 
dicator of the magnetic behavior of the control samples. 
The substitutional fraction for the control samples for 
regions-A and C shown in Fig. [2] is larger than the sam- 
ples with a Pb-surfactant. It is likely that these Mn s „j 
are also substitutional-like sites within the precipitates. 
Support for this idea is found in the results of MBE grown 
Sii-aMiXj, alloys prepared by co-deposition of Mn and Si, 
where the alloy undergoes spinodal decomposition into 
nanocolumnsi^ For 1.5% Mn, RBS shows a comparable 
x S ub = 0.7 at Tsi = 200 °C. We have recently repro- 
duced these growth conditions and also observed the for- 
mation of nanocolumns. However, despite the large x su b 
reported in Ref. [H, these nanocolumns are paramagnetic 
with a small moment^ 

To understand how Pb can influence the formation of 
Mngi we performed DFT calculations. In our study, we 
find that Pb increases the formation energy of Mnsi by 
direct substitution into a surface Si site. However, the 
Pb surfactant tends to weaken the top layer bonding of 
Si because of the overly crowded bonding to the top Si 
atom. Our DFT calculations show that this enhances 
vacancy formation in the top Si layers, and thereby 
creates a channel for the Mn to incorporate substitu- 
tionally. Various theoretical and experimental studies 
show that different amounts of Pb on Si can form differ- 
ent surface reconstructions 5 -^—. For one full monolayer 



coverage of Pb, we picked a 2 x 1 surface reconstruc- 
tion as both confirmed by our in situ RHEED obser- 
vations and by scanning tunneling microscopy measure- 
ments and calculations 5 -^—. We performed DFT calcu- 
lations using the VASP code within the generalized gra- 
dient approximatio n 54 ! 55 . We chose the Projected Aug- 
mented Wave PBE potential as our pseudo potentia l 56 ' 57 . 
Our simulation cell is a 144 Si atoms slab with 32 H atoms 
passivating the dangling bonds of the bottom surface as 
shown in Fig. [SJa). One full layer of Pb atoms were 
introduced into the calculation as the surfactant. The 
thickness of the vacuum layer is about 18 A. We used a 
plane wave cutoff energy of 245 eV and a 2 x 2 x 1 k- 
point mesh for Brillouin zone sampling. A total energy 
minimization was performed by relaxing atomic positions 
until the forces converged to less than 0.01 eV/A. The 
calculated Si lattice constant is 5.469 A, which is in good 
agreement with the experimental value of 5.43 A. 




(a) (b) 



FIG. 8. (Color online) (a) Schematic illustration of the sim- 
ulation supercell. Middle (grey): Si. Bottom small (white): 

H. (b) Schematic illustration of the surface reconstruction of 
one monolayer of Pb on Si thin film. Top big (black): Pb. 
Middle: Si. Bottom small (white): H. 

We define the formation energy of a vacancy as 
E '(formation) = E(vac) — E (reference) + fi(Si), where 
E (reference) is the total energy of the Si slab, E(vac) 
is the total energy of the same cell with one Si atom 
removed and fJ.(Si) is the chemical potential of the Si, 
which is the energy per atom of a bulk calculation of a 
system with 512 Si atoms. The chemical potential of the 
Si is calculated to be —5.43 eV. We calculated the forma- 
tion energy of vacancy of the atom on the surface layer 
in a pure Si thin film by removing a Si atom from the 
top dimer layer of the 2x1 reconstructed surface. We 
found that the vacancy formation energy in this case is 

I. 18 eV. Next, we introduced Pb atoms as surfactants 
into the system, as shown in Fig. [8] (b). We removed one 
Si atom from the top Si layer and calculated a vacancy 
formation energy of 0.26 eV, which is about 0.92 eV lower 
than the vacancy formation energy for a pure Si surface 
at the same site. This is mainly due to the Pb surfac- 
tant. The Pb atoms largely disturb the top Si surface 
layer and lower the Si vacancy formation energy. On the 
other hand, the interstitial diffusion is unlikely to change 
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much due to the surfactant. 



V. CONCLUSION 

Our results demonstrate the control of the incorpora- 
tion of Mn 5-doped layers into Si(OOl) via a Pb surfac- 
tant. We find that the resulting structure is sensitive to 
the growth conditions of the Si capping layer. Structural 
and magnctomctry studies identify three distinct regions 
of the magnetic phase diagram, where two of them are 
ferromagnetic. Samples in rcgion-A of Fig. [1] have a Tc 
above 400 K and contain nanocrystalline MnSi precip- 
itates with a B2-like crystal structure. This high tem- 
perature ferromagnetism extends into region-C, where 
a second magnetic transition at Tc = 170 K exists 
due to MnSi nano-disks with a B20-like crystal struc- 
ture. This result complements previous work on MnSii.7 
precipitates^ ' 8 ' 9 ' 11 ' 12 and provides examples of different 
nanocrystalline precipitates in Si. More work needs to be 
done with MnSii.7 and MnSi-B20 to determine the role 
that magnetic defects play in these structures.— 

In this paper, we discuss a mechanism for an intrigu- 
ing surfactant enhanced substitutional doping phenom- 
ena. In the two ferromagnetic phases, RBS shows that 
an increase in m r is possibly correlated with an increase 
in substitutional-likc Mn. Wc interpret this increase as 



an increase in the concentration of substitutional Mn in 
the Si matrix, which may be the origin of Tq > 400 K 
phase. DFT calculations show that a Pb layer on Si re- 
duces the vacancy formation energy by 0.92 eV which 
opens the pathway for substitutional doping of Mn into 
Si. This provides a new opportunity for doping in other 
DMS materials. 
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